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Abstract 

We introduce Eulerian maps with blocked edges as a general way to imple- 
ment statistical matter models on random maps by a modification of intrinsic 
distances. We show how to code these dressed maps by means of mobiles, i.e. 
decorated trees with labeled vertices, leading to a closed system of recursion 
relations for their generating functions. We discuss particular solvable cases in 
detail, as well as various applications of our method to several statistical sys- 
tems such as spanning trees on quadrangulations, mutually excluding particles 
on Eulerian triangulations or the Ising model on quadrangulations. 



1. Introduction 

Enumeration of maps, i.e. connected graphs embedded in a surface, is a funda- 
mental issue that has drawn the attention of mathematicians and physicists all over the 
years ever since the seminal series of papers of Tutte on the subject [Q. In physical 
applications, the subject was extended so as to include statistical models defined on 
the maps, by equipping them with "matter" configurations, say of particles, dimers, 
spins or trees. These dressed maps were introduced as discrete realizations of random 
surfaces in the general framework of two-dimensional quantum gravity (2DQG) for 
which many results have been obtained so far (see for a review). Several techniques 
have been used to address this problem of counting maps, with their possible mat- 
ter decorations generally weighted by Boltzmann factors. Beside the original recursive 
decomposition technique of Tutte, mainly applied to maps without matter, the most 
powerful approach was undoubtedly that based on matrix integrals 0]. Many results 
obtained that way display a simplicity contrasting with the high sophistication of the 
method. In particular, most planar enumeration results involve algebraic generating 
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functions. This led combinatorists to look for alternative and more constructive ap- 
proaches that could explain this simplicity Quite recently, a general technique 
was developed that relies on bijections between planar maps with possible decorations 
and families of decorated trees . The bijections are of two different types. A first class 
of bijections well adapted to maps with prescribed vertex valences, consists of a 
cutting algorithm of the map into a blossom tree that carries the minimal information 



needed to recover the map. A second class of bijections |1TQ,11| is well adapted to the 



dual version of the former, namely maps with prescribed face valences. It makes use of 
the intrinsic geodesic distance on the map to code it by a so-called labeled mobile, i.e. 
a tree decorated by integer labels recording this distance in a way that allows to recon- 
struct the original map. Both methods were used to enumerate families of bicolored 
planar maps, or equivalently their dual Eulerian maps, in the absence of matter, all dis- 
playing the standard "square-root singularity" characteristic of pure gravity (if we omit 
non-generic, non-physical multicritical points). Still, one may in some cases reinterpret 
the color as a matter degree of freedom via suitable weightings. This gives access for 
instance to the solution of the Ising model on tetravalent maps [0, leading to a new 
singular behavior of the generating function at the price of some analytic continuation 
to reach the transition point. Similarly, this allows to treat the case of hard particles 
on arbitrary (non bicolored) maps, but with no new singular behavior in this case . 



In Ref . , it was found how to extend the bijection with blossom trees to the case 
of bicolored maps equipped with interacting particles, allowing in particular to recover in 
a purely combinatorial way the crystallization transition of hard particles. Considering 
more involved particle exclusion rules leads to most of the possible universality classes 
of rational conformal theories (RCFT) with central charge c < 1 coupled to 2DQG ||15|| . 

The aim of this paper is to extend the mobile formalism to this same family of 
bicolored (or dually Eulerian) maps with mutually excluding particles. As we shall see 
however, our method proves more general and applies to other classes of matter such 
as spanning trees or Ising spins, here without analytic continuation, thus providing 
a unified framework for all these models. More precisely, as already mentioned, the 
original mobile coding of Eulerian maps uses the geodesic distance from a given origin 
vertex on the map. A way to incorporate matter is precisely to modify this distance 
by introducing blocked edges, i.e. edges excluded from geodesic paths. All the matter 
systems above can indeed be reformulated in terms of configurations of blocked edges 
on Eulerian maps. The mobile construction is then easily modified so as to account for 
the presence of the blocked edges and this eventually leads to an efficient coding of the 
maps at hand in the form of labeled trees with marked edges. Using the tree structure, 
one can derive recursive equations for the corresponding generating functions, leading 
to a combinatorial solution of the original problems. 

The paper is organized as follows. Section 2 is devoted to the general theory of map 
coding via mobiles. We first define precisely in section 2.1 the notion of Eulerian maps 
with blocked edges and show in section 2.2 how to code them by labeled mobiles that 
we characterize in detail. The procedure leading back from the mobile to the original 
map is described in section 2.3. In section 2.4, we further use the characterization of 
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mobiles to obtain a system of recursive equations for their generating functions, also 
easily translated into map generating functions. 

Sections 3 and 4 describe two large classes of applications. In section 3, a first 
series of direct applications consists of models where the blocked edges themselves form 
the matter degrees of freedom. Section 3.1 deals with the statistics of Eulerian maps 
with bi- and tetra-valent faces and with blocked edges, which can also be viewed as 
decorated quadrangulations whose unoriented edges can be blocked in both directions 
independently, like maps with possible one-way roads and/or fully blocking roadworks. 
The simpler case of quadrangulations with fully blocked edges only is treated in section 
3.2 where a detailed analysis of the critical singularities is presented. This allows in 
particular to solve the problem of quadrangulations with a spanning tree. In appendix 
A, we provide for both solutions a simple combinatorial interpretation in terms of maps 
with vertices of even valence and without matter. Appendix B deals with the case of 
maximally blocked Eulerian maps, i.e. maps with a maximal number of blocked edges 
but which remain connected in the sense that every vertex can still be reached via some 
path from the origin. 

Section 4 is devoted to applications to maps with particles subject to exclusion 
rules. We first show in section 4.1 how to formulate these problems in terms of Eulerian 
maps with blocked edges by associating to any unconstrained configuration of particles 
on a map a set of configurations of blocked edges on the same map with possible negative 
weights, and whose total contribution vanishes as soon as the particle exclusion rule is 
violated. This allows for the derivation in section 4.2 of recursive relations for the 
generating functions of the models. Two particular examples are discussed: Eulerian 
triangulations with hard particles in section 4.3 and, in section 4.4, the Ising model on 
quadrangulations via a straightforward bijection with hard particles. Wc gather a few 
concluding remarks in section 5. 



2. Coding maps with mobiles 

Eulerian maps are planar embedded graphs (considered up to continuous deforma- 
tion), which may be drawn in a single loop without lifting the pen. As each vertex 
is of even valence, the edges of these maps are naturally oriented, with incoming and 
outgoing edges alternating around each vertex. We may then imagine a walker moving 
along these oriented edges as on one-way streets throughout the map. We wish here 
to introduce the notion of Eulerian map with blocked edges, by which we simply mean 
that some of the edges of the map may be forbidden to the walker, like streets blocked 
by road works or so. The precise position of the blocked edges will be in general arbi- 
trary up to one important restriction: we demand that every vertex of the map may be 
reached from a chosen origin vertex despite the blockings. As we shall see later, many 
statistical models on random maps may be described in terms of such Eulerian maps 
with blocked edges. This includes for instance spanning trees, hard particles, or even 
Ising spins. 
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2.1. Eulerian maps with blocked edges 

Let us now come to more precise definitions of Eulerian maps witfi blocked edges. 
An Eulerian map is equivalently defined as a planar map wliose faces are colored in 
black and white so that no two adjacent faces be of the same color. In other words, its 
dual mapEI is a bipartite planar map, namely a map whose vertices may be partitioned 
into two subsets such that no two vertices of the same subset are adjacent. The coloring 
of the faces of an Eulerian map induces a natural orientation of its edges by simply 
demanding for instance that this orientation be clockwise around black faces. Here we 
shall consider so-called pointed maps, i.e. maps with a distinguished vertex, hereafter 
referred to as the origin. 




Fig. 1: An example (a) of pointed Eulerian map with blocked edges (thick 
lines). The blockings satisfy the global connectivity constraint that every 
vertex can be reached from the origin by a path consisting of non-blocked 
edges taken in their canonical orientation (indicated by an arrow). On 
the dual bipartite map (b) , the duals of blocked edges necessarily form a 
forest, i.e. a graph without cycles and, moreover, are such that each face 
can be reached from the external face (dual to the origin) without crossing 
the forest edges and by going clockwise (resp. anticlockwise) around black 
(resp. white) vertices. 

By map with blocked edges, we simply mean that we restrict the set of paths allowed 
on the map by deciding that some of the edges of the map are "blocked" , while paths are 

^ Recall that the dual of a planar map is the planar map obtained by replacing each face 
by a vertex at its center and each edge by a "perpendicular" edge connecting the centers of 
the adjacent faces. The vertices of the original map become the faces of the dual map. 
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only allowed to pass through un-blocked edges and must respect the orientation on these 
edges. As mentioned above, we furthermore impose a global connectivity constraint that 
each vertex be attainable from the origin by some allowed path (see figure l-(a) for an 
example). From this requirement, we immediately see that, on the dual bipartite map, 
the set of edges dual to blocked edges (also referred to as blocked edges for convenience) 
form a forest, i.e. contains no cycle. This forest is however not arbitrary in general 
as it must be such that any face of the dual map can be reached from the origin face 
(dual of the origin vertex of the Eulerian map) without crossing blocked edges and by 
going clockwise (resp. anticlockwise) around black (resp. white) vertices (see figure 
l-(b)). The connectivity constraint is crucial for the construction below to be valid. 
In the following, we shall be led to introduce further restrictions on the blocked edge 
configurations, whose consequence on our construction will be transparent. 

2.2. Well-labeled mobile construction 




(i) (ii) (iii) 

Fig. 2: The rules for the construction of a well-labeled mobile (see text). 
A fiagged edge can be marked (thick line) or not (thin line) according to 
whether the dual edge was blocked (i) or not (ii). 

The connectivity constraint above allows to assign to each vertex of the map its 
(finite) distance from the origin defined as the minimal number of steps of a path going 
from the origin to that vertex via un-blocked edges only and respecting the orientation. 
Such a path will be referred to as a geodesic path. 

We can encode the configuration of the map and its blocked edges via a so-called 
well-labeled mobile as follows: 

1. We first label all vertices of the map by their distance from the origin. 

2. We then add a white (resp. black) vertex at the center of each white (resp. black) 
face. 

3. We then consider each edge of the map and apply the following construction: 

(i) If the selected edge is a blocked edge, we draw its dual edge by connecting 
the black and white vertices at the center of its adjacent faces and complete 



5 



Fig. 3: Construction of tiie well-labeled mobile associated with the Eule- 
rian map with blocked edges (a) of figure 1. Applying the rules of figure 
2 for each edge of the map (b) results in the desired mobile (c). On this 
mobile, flagged edges can be marked (thick lines) or not. 



it with a flag on each side. We put on each flag the label of the closest vertex 

(see figure 2-(i)). We finally mark this dual edge. 
If the selected edge is not blocked, it points from a vertex at distance n from the 
origin to a vertex at distance m, with necessarily m < n + 1. 

(ii) If m < n, we repeat the above construction and draw the dual edge by con- 
necting the black and white vertices at the center of its adjacent faces and 
complete it with a flag on each side, with again the same label as that of the 
closest vertex (see flgure 2-(ii)). The dual edge is however un- marked in this 
case. 

(iii) If m = n -|- 1, we draw an edge from the white vertex at the center of the 
adjacent white face to the vertex labeled m = n + 1 (see flgure 2- (iii)). 

The well-labeled mobile is simply deflned as the graph obtained from the collection of 
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all the black and white vertices, all the added marked (case (i)) and un- marked (case 
(ii)) flagged edges with their labeled flags, and all the added edges from white vertices 
(case (iii)) together with their labeled endpoints (see flgure 3 for an example). For 
convenience, we shall refer to marked flagged edges, un-marked flagged edges and edges 
between white vertices and labeled vertices as edges of type (i) , (ii) and (iii) respectively. 
Note that, by construction, all labels on labeled vertices are positive integers (since the 
origin vertex cannot be in the mobile), while all labels on flags are non- negative integers. 
The above construction extends that of Ref. [l^] that associates well-labeled mobiles to 
Eulerian maps without blocked edges by use of (ii) and (iii) only. 

Let us now come to the characterization of the well-labeled mobiles obtained from 
this construction. The flrst and most important characterizing property of a well-labeled 
mobile is that it is a tree, i.e. a connected graph without cycles. The proof is similar 



to that of Ref. [|Tl|]. First we note that edges belonging to a geodesic path yield only 
un-marked edges of type (iii) in the mobile, as the distance increases by one at each 
step. This guarantees that geodesic paths on the map never cross the flagged edges 
of the mobile, and may cross the mobile only at its labeled vertices. We then assume 
by contradiction the existence of a cycle on the mobile, and call interior the region 
delimited by this cycle that does not contain the origin. We then consider the smallest 
label n among the labels of the (labeled) vertices on the cycle and the labels of the flags 
along the cycle that lie in the interior region. If n is attained for a labeled vertex on 
the cycle, we immediately deduce from the construction (iii) the existence of a vertex at 
distance n — 1 from the origin in the interior. Any oriented geodesic path from the origin 
to that vertex must intersect the cycle at a labeled vertex with label strictly smaller 
than n. This is a contradiction. If n is not attained at a labeled vertex, but only at 
a flag, we deduce from the construction (i) or (ii) the existence of vertex at distance n 
from the origin in the interior, hence a geodesic path to that vertex intersects the cycle 
at a labeled vertex with label < n. This is again a contradiction. We deduce that the 
mobile has no cycle, hence is a forest made of c connected components, where c is also 
the difference between its number of nodes and its number of edges. By construction, 
this number of edges is simply the number E of edges of the original map, while the 
number of nodes is equal to the number of faces F of the original map plus the number 
of labeled vertices. The latter is at most y — 1 if we denote by V the total number 
of vertices of the original map, as the origin cannot belong to the mobile (it remains 
isolated by construction hence is removed) . AsF + V — E = 2 from the Euler relation 
for a planar map, we deduce that c<F + {V — 1) — E = 1, hence necessarily c = 1. 
The mobile is thus a tree that moreover contains all the vertices of the original map but 
the origin. 

The second characterizing property of well-labeled mobiles concerns the environ- 
ment of white and black vertices, i.e. the conflguration of labels around them. Recall 
that, by corner at a vertex is meant any sector delimited by two consecutive incident 
edges around that vertex. In the case of a white vertex, we flrst note that, by construc- 
tion, the number of incident edges is nothing but the valence of the associated white 
face on the original map. We have the following rules for the cyclic evolution of the 
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Fig. 4: Constraints on tlie labels around a white vertex of a mobile (see 
text). 

labels around each white vertex (see figure 4): 

- at the crossing of each un-marked fiagged edge, the value of the flag must increase 
weakly clockwise. There is no such restriction for marked flagged edges. 

- going through each corner clockwise around a white vertex, the label decreases by 
one if the first edge of the corner is of type (iii), and remains constant if the first 
edge is of type (i) or (ii). 




Fig. 5: Constraints on the labels around a black vertex of a mobile (see 
text) 

As for black vertices, their incident edges are of type (i) and (ii) only and we have the 
following rules for the cyclic evolution of the labels around each black vertex (see figure 

8 



5): 

- at the crossing of each un-marked flagged edge, the value of the flag must decrease 
weakly clockwise. There is no such restriction for marked flagged edges. 

- going through each corner clockwise around a black vertex, the label increases 
weakly. 




Fig. 6: An alternative representation using spurious dangling edges of the 
constraints of figure 5 for the evolution of labels at a corner around a black 
vertex. 

Note that the first (crossing) rules for white and black vertices are actually redundant 
with each other, but we mention both of them to emphasize the symmetry. Moreover, 
as displayed in figure 6, the second (corner) rule for black vertices can be made similar 
to that of white vertices by dressing each corner around a black vertex with first label n 
and second label m > n clockwise by m — n spurious dangling (wiggled) edges attached 
to spurious vertices labeled 7i + l,7i + 2,...,m. With this dressing, the black corner 
rule becomes: 

- going through each corner clockwise around a black vertex, the label increases by 
one if the last edge of the corner is a dangling edge and remains constant if the last 
edge is of type (i) or (ii). 
In this latter formulation, the total number of edges around a black vertex, including 
the spurious ones, is nothing but the valence of the associated black face on the original 
map. 

The above properties for the evolution of labels, together with the fact that these 
labels are positive integers on vertices and non-negative integers on fiags, with at least 
one vertex labeled 1 or one fiag labeled 0, form a complete characterization of the well- 
labeled mobiles. In the following we will be led to relax the above constraints on labels. 
Mobiles for which we relax the constraint of existence of a label 1 or but keep the 
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positivity constraint will be referred to as positive mobiles. Finally, mobiles for which 
we relax both constraints will simply be referred to as unrestricted mobiles. 



2.3. Inverse construction 

The well-labeled mobiles constructed above form compact encodings of the asso- 
ciated Eulerian maps with blocked edges. Indeed, to recover the original map and 
its blocked edges from the corresponding mobile, we simply apply the construction of 
Ref.lTTI] as follows. First, by labeled corner we simply mean a corner at a labeled vertex 



of the mobile, with the corresponding label value. We then move clockwise around the 
mobile and record the cyclic sequence of successive labels encountered at flags and/or 
labeled corners (spurious edges and vertices displayed in flgure 6 are not considered). 
The passage from a label to the next is always via a white or a black vertex, hence, 
from the rules obeyed by labels around those vertices, we immediately deduce that the 
sequence of labels decreases by one after each labeled corner and increases weakly after 
each flag (in particular, it remains constant when passing along a white vertex). Oth- 
erwise stated, the sequence of labels around the mobile has the "ratchet" property that 
it can have arbitrary increasing steps but decreasing steps of at most 1 which take place 
just after labeled corners. 

This property immediately ensures the existence of successors deflned as follows: 
the successor of a labeled corner with label n > 1 is the flrst labeled corner encountered 
clockwise with label n — 1. Similarly, the successor of a flag with label n > is the 
flrst labeled corner encountered clockwise with label n. Note that the ratchet property 
of the label sequence also guarantees that if a corner or flag /' lies between a corner or 
flag I and its successor s(/), the successor s(l') is necessarily attained at or before s{l). 
We may therefore connect each labeled corner with label n> 1 and each flag with label 
n > to its successor by a new edge, called a chord, in such a way that all these chords 
do not cross, while all corners labeled 1 and all flags labeled remain adjacent to the 
exterior face. We flnally add an extra origin vertex inside the exterior face and connect 
all corners labeled 1 and all flags labeled to that vertex by non-crossing chords. The 
original Eulerian map is made of the origin, the labeled vertices and the chords. Note 
that some chords (those starting from flags) have to be glued at the level of the double 
flags so as to form the desired edges between the (origin or labeled) vertices of the map 
(see flgure 7 for an example). Any edge of the Eulerian map is of one of the following 
three types (see flgure 7): 

(i) either it is obtained by gluing two chords at the level of the double flag of a marked 
edge. We then make it a blocked edge. 

(ii) or it is obtained by gluing two chords at the level of the double flag of an un-marked 
edge. This edge is then left un-blocked and we decide to orient it so that the black 
endpoint of the flagged edge lies on its right. 

(iii) or it is obtained from a single chord connecting a labeled vertex to its successor. 
This edge is then left un-blocked and we decide to orient it backwards from the 
successor. 

The proof that the above construction and that of the previous section are inverse of 
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Fig. 7: Inverse construction: starting from a well-labeled mobile (a) , each 
labeled corner or flag is connected to its successor (b). The graph formed 
by these chords (c) is the desired Eulerian map (d), where the blocked 
edges are obtained from the chords originating from marked flagged edges. 



one another is similar to that of Ref. |Tl|] in the absence of blocked edges and will not 
be reproduced here. The orientation of the edges simply reproduces that induced by 
the black and white faces on the Eulerian map. Finally, the labels on the vertices 
simply record the distance of these vertices form the origin by oriented paths avoiding 
blocked edges. Indeed, form our choice of orientation above, the sequence of labels along 
any oriented path from the origin on the map either decreases weakly (case (ii) above) 
or increases by one (case (iii) above) at each step. The label of a vertex is therefore 
necessarily smaller than its distance from the origin. In other words, the distance from 
the origin to a vertex is larger or equal to its label, and it must be exactly equal as the 
sequence of consecutive successors of that vertex provides, when considered backwards, 
an oriented path from the origin to that vertex whose length is exactly equal to its label. 
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2.4- Generating series 



The main interest of the above coding by mobiles concerns the enumeration of 
famihes of Eulerian maps with blocked edges. Indeed, such enumeration simply amounts 
to counting appropriate mobiles and the underlying tree structure of these mobiles 
translates into simple recursive equations for the associated generating functions. 

In the following, we shall consider only mobiles with at least one labeled vertex. 
In other words, we exclude the degenerate case of mobiles made only of flagged edges, 
necessarily all with equal labels. Going back to maps, this amounts to excluding the 
quite trivial maps made of a single vertex only, connected by a configuration of nested 
"petals" . The enumeration of mobiles is then made easier by considering rooted mobiles, 
i.e. mobiles with a distinguished corner around a labeled vertex. We may then view the 
mobile as hanging from an extra root edge added in this corner. 

Let us now consider the quite general problem of enumerating maps with prescribed 
face valences and with a fixed number of blocked edges. At this stage, we make for sim- 
plicity no restriction on the configuration of blocked edges around a face other than 
the global connectivity constraint. Some extra restrictions will be introduced in some 
applications below, whose incorporation will lead to straightforward and harmless mod- 
ifications. As usual, it is simpler to work within grand canonical ensembles of maps 
in which the numbers of faces of given valences may vary and are controlled by weight 
factors. More precisely, we consider the set of maps with faces of arbitrary valences up 
to some upper bound p, counted with weight gk per /c-valent white face, gk per /c-valent 
black one {k — 1,2, . . . ,p). We also let the number of blocked edges vary, with a weight 
y per blocked edge. This amounts to consider mobiles with white and black vertices 
of valence k < p (including the spurious dangling edges for black vertices), weighted 
respectively by gk and gk, and with an arbitrary number of marked edges weighted by 



We denote by Rn the generating function of such weighted mobiles rooted at a 
corner labeled n. At this stage, we need not specify whether the labels obey positivity 
conditions (case of positive mobiles) or not (case of unrestricted mobiles). Indeed, the 
following reasoning holds in both cases, the difference coming only from boundary con- 
ditions that we shall discuss later. We shall not consider the case of well-labeled mobile 
(rooted at a corner labeled n) as their generating function is directly expressed in 
terms of that for positive mobiles via Gn = Rn — Rn-i- 

Let us derive the set of recursion relations for Rn inherited from the mobile struc- 
ture. We have the first relation 



where Ln is the generating function of mobiles rooted at a univalent vertex labeled n. 
This is readily seen by decomposing the root vertex according to its arbitrary valence 
k (see figure 8). 




(2.1) 
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Fig. 8: Recursive decomposition of a mobile rooted at a labeled corner. 
When the root vertex has valence k, the decomposition yields k mobiles 
rooted at a univalent labeled vertex. In terms of generating function, this 
translates into eq.(2.1). 



R, 



B 
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Fig. 9: The three possible types of descending subtrees at a white vertex 
of a mobile according to the type of their root edge, together with the 
associated generating functions. 



The generating function L„ is then obtained by inspecting the environment of the 
white vertex attached to the univalent root vertex. Its descendent subtrees come into 
three species (see figure 9): 

(i) sub-mobiles rooted at a corner around a labeled vertex, with generating function 
Ri if the root label is i, 

(ii) sub-mobiles rooted at some un-marked flagged edge incident to a black vertex, with 
generating function Bij if i (resp. j) is the label on the right (resp. left), viewed 
with the white vertex on top, 

(iii) sub-mobiles rooted at some marked flagged edge incident to a black vertex, with 
generating function Bi j if again i (resp. j) is the label on the right (resp. left). 

To follow the evolution of labels clockwise around the inspected white vertex, we in- 
troduce a formal (transfer) operator Qo acting on basis vectors \i) indexed by the 



13 




j<i 

Fig. 10: The three possible types of descending subtrees at a black ver- 
tex of a mobile according to the type of their root edge, together with 
the associated generating functions. It proves convenient to use the rep- 
resentation of fig. 6 of mobiles with spurious dangling edges (with a trivial 
generating function 1) which emphasizes the similarity between black and 
white vertices. 

encountered labels. The action of this operator reads: 

Qo\^) = R^\^-l) + , (2.2) 

where the ranges of summation incorporate the white crossing rules (obeyed by labels 
at the crossing of flagged edges) displayed in fig. 4. Returning to the computation of Ln, 
we note that, as the first label encountered clockwise is n — 1, the succession of labels 
around the white vertex must lead from n — 1 back to n in /c — 1 steps if the white vertex 
is /c-valent. This leads to 

Lr, = Y,gk{n\Cll-^\n-l) , (2.3) 

k 

where {i\ denotes a vector of the dual basis defined by = 5i^j. 

Likewise, inspecting the environment of a black vertex incident to a flagged root 
edge, we get 

k 

and 

B^,J=yJ29k{m':-'\^) , (2.5) 

k 

where the transfer operator Q, around a black vertex now reads 

Q.|z) = \^ + l) + J2 W,,j\j) + J2 W^,J\J)■ (2.6) 
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Here the three terms account for the three species of sub-mobiles one may encounter 
around a black vertex (see figure 10): 

(i) dangling ends with weight 1, 

(ii) sub-mobiles rooted at some un-marked flagged edge incident to a white vertex, with 
generating function Wij if i (resp. j) is the label on the right (resp. left), viewed 
with the black vertex on top, 

(iii) sub-mobiles rooted at some marked flagged edge incident to a white vertex, with 
generating function Wij if again i (resp. j) is the label on the right (resp. left). 

The ranges of summation in eq.( p.6| ) now incorporate the black crossing rules (obeyed 
by labels at the crossing of flagged edges) displayed in flg.5. The generating functions 
Wij and Wij are flnally obtained from: 

W.,j=J29kmo-'\^) . J<^ (2.7) 

k 

and 

m,j=yJ29km':-'\t) . (2.8) 
k 

Let us now discuss the distinction between positive and unrestricted mobiles. In the 
case of positive mobiles, the equations (|2.1|) - (p.8| ) above are only valid for the allowed 
range of labels, i.e. n > for and and i,j > for Bij, Bij, Wij and Wij. 



Moreover, in eqs. ( p.2|) and (|2.6| ), we must implicitly assume that Ri = for i <0 and 
Bi j = Bi j = as well as Wij — Wij = if i < or j < 0. With these conditions, 
the set of equations ( p.l| )-( PTBD determine completely all generating functions as formal 
power series of the g^s and the ^^'s. 

The case of unrestricted labels is much simpler as (|2.iP -( pT8| ) are valid for all values 
of the indices and moreover, enjoy a manifest translation invariance property. More 
precisely, we may write 

R = Rn 

> 

(2.9) 



R = 


Rn 




Bi = 


Bn, 


n+e 


Be = 


Bn, 


n+e 


We = 


Wn 


,n+e 


We = 




,n+£ 



as all the above quantities on the right hand side are independent of n. 

Going back to maps, one can easily check that R is the generating function for 
Eulerian maps with blocked edges, with an origin vertex and with a distinguished edge 
which is not blocked and moreover points from a vertex at distance m from the origin to 
one at distance m-fl, for some m. Indeed, we have clearly a bijection between the set of 
well-labeled mobiles rooted at a corner with arbitrary label, and the set of unrestricted 
mobiles rooted at a corner with prescribed label (say 0), whose generating function is 
R. The bijection simply consists in shifting all labels of the well-labeled mobile by a 
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constant value, such that the root label becomes 0. Finally, the corners labeled m + 1 on 
a well-labeled mobile are in correspondence with non-blocked edges of type m — > m + 1 
on the corresponding map. 

One may as well distinguish an edge on the map which is not blocked, but now 
points from a label at distance m from the origin to a label at distance m — i {£ > 
0) for some arbitrary m. The corresponding generating function is easily seen to be 
BgW-e- Finally, one may also distinguish an edge which is blocked and connects a 
vertex at distance m from the origin to a label at distance m — i for some arbitrary 
m. The corresponding generating function reads {l/y){B£W-£ + B-eWe) for £ 7^ and 
{l/y)BoWo for i = 0. To summarize, the generating function of pointed Eulerian maps 
with blocked edges and with a distinguished non-blocked edge is 



R + J^BeW-i (2.10) 



£>0 



while that of pointed Eulerian maps with blocked edges and with a distinguished blocked 
edge is 

. (2.11) 

2/ V 

The functions i?. Be, Wi, Bi and Wi are the translation invariant solutions of 
the equations (|2.1| )- (|2.8| ), which can be rewritten in a slightly more compact form by 
introducing a formal variable z and defining the Laurent series Qdz) and Q»{z) via 

Qoiz) = - + J2Biz' + J2B£/ 

' (2.12) 



Q.{z) = z -fJ^W.^ + ^W-,. 



that are translation invariant versions of eqs.( |2.2| ) and (|2.6|) . The recursive equations 
read: 



k 



B, = Y,~9kQ.{zf-\^, ^>0 

k 

Bi = y ^gkQ»{zf~^\^, (2.13) 

k 

k 

k 

where stands for the extraction of the coefficient of z^ in the corresponding Laurent 
series. 
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Eliminating L, all S's and W^s, we arrive at the closed system of equations 



k k 

Q.[z)^z + Y.9k [Ql~\z)]_ + yY,9kQ':-\z) 

k k 



R=l + J29k 

k 



R 



(2.14) 



where the notations [-j-i-, [■]_ and [-jo stand respectively for the extraction of the non- 
negative power part, the non-positive power part, and the constant term in the Laurent 
series. In practice, this system may be solved in two steps: the first two lines determine 
all coefficients of the Laurent series as functions of R and formal power series of the 
(7's and ^'s. Substituting these in the third line yields a self-consistent equation for R, 
which itself is solved order by order in the (7's and ^'s. 

When y = we recover the usual equations for Eulerian maps (without blocked 
edges), as obtained from the mobile construction of Ref. [11 
thogonal polynomial solution of the two-matrix model. 



or by the so-called or- 



It is interesting to note that 
the restriction of Laurent series to their positive or negative parts simply reflects the 
constraints on flag labels for un- marked edges (crossing rules). For marked edges (cor- 
responding to the y term) there is no such restriction. 

Finally, we have the duality property that, whenever gk = Qk for all k, Qo{z) is 
equal to (5,(— ), while R can be alternatively recovered via: 



(2.15) 



3. Direct applications 

In this section, we show for illustration how to use the above formalism to study 
the statistics of blocked edges in a number of simple cases in which the blocked edges 
themselves form the matter degrees of freedom on the maps. 

3.1. Blocked edges on quadrangulations viewed as Eulerian maps 

As a first direct application of the results of section 2, let us consider the simple 
case of quadrangulations viewed as Eulerian maps as follows: consider Eulerian maps 
where all white faces are tetravalent and all black faces are bivalent. We may squeeze 
all the black faces into single edges, which results in a (non necessarily Eulerian) map 
with tetravalent faces only, i.e. a quadrangulation. After this transformation, every 
edge replaces two edges of the original map, with opposite orientations. Hence, in 



17 



Fig. 11: An example of pointed Eulerian map (left) witli blocked edges 
(thick lines) satisfying the global connectivity constraint, where all black 
faces are bivalent and all white faces tetravalent. Each black face can 
be squeezed into an edge, resulting in a quadrangulation (right) whose 
edges can be either non-blocked (thick lines) , blocked in one direction only 
(thick lines with one-way signs) or blocked in both directions (thin lines 
with roadwork signs). 

the absence of blockings, each edge of the quadrangulation can now be taken in both 
directions. The introduction of blockings now corresponds to forbidding either one 
direction (one way road) or both, keeping the global connectivity constraint that every 
vertex be reachable from the origin of the map (see figure 11 for an example). Taking 
eq.( p.l4|) with — gSk,4 and gk = Sk,2 (as there are exactly twice as many black 
bivalent faces as white tetravalent ones, we need not introduce an extra weight factor 
for the black faces), we get the system: 



Qo{z) 

R = l+g 



R 

z 
z 



- + [Q.{z)]+ + yQ,{z) 

+ 9 [QUz)]_+9yQUz) 



(3.1) 



First we note that, by parity, both Qo{z) and Q,{z) contain only odd powers of z. From 
the second line, we then immediately deduce that 



[Q.{z)]^ = z + gy [QUz)]^ 
and, upon substituting in the first line 



R 

Qo{z) = - + {l + y)z + gy{l + y)Qo{zf . 
z 



(3.2) 



(3.3) 
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Upon setting x = gy{l + y){R/z + (1 + y)z)'^ and Qoiz) = {R/z + (1 + y)z)q{x), this 
equation becomes 

q{x) = 1 + xq{x)^ 

whose formal solution is 

(3?i)! 



q[x) 



Writing the last line of eg. (PH]) as 



R = l+g 



1 



{3n)\ 



E 

n>0 



R 



n!(2n + 1)! 



X 



Qo{z) {l + y)z]- 

z z 



R 



J 



r^/" (2/(1 + 2/)) 



n-l 



R 



z 



+ (1 + 2/)^ 



2n+l 



^ n!(2n + 1)! 

n>l ^ 

we end up with a single equation for the generating function R 



R 



r=i+j: 



(3n)! 



^ n!n!(n + 1)! 



(3.4) 
(3.5) 



(3.6) 



(3.7) 



Rather than directly analyzing this formula, we will derive and study in the next section 
a very similar and slightly simpler formula (see eq.(3.10) below). Their connection is 
discussed in appendix A below. 



3.2. Blocked edges on quadrangulations and spanning trees on the dual 



Fig. 12: A quadrangulation (left) with edges (thick lines) blocked in both 
directions, and satisfying the global connectivity constraint (irrespectively 
of the choice of the origin vertex). On the dual tetravalent map (right), 
duals of blocked edges (thick lines) form a forest. 
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As there is no canonical orientation of the edges in quadrangulations, a more nat- 
ural definition of blocked edges consists in either allowing for both directions on a given 
edge or forbidding both simultaneously. In the Eulerian map equivalent formulation of 
previous section, this simply means that blockings always take place for both directions 
simultaneously. In this case, for a fixed map, the configurations of blocked edges sat- 
isfying the global connectivity constraint turn out to be independent of the choice of 
the origin. Rephrased in the dual language, this constraint simply amounts to demand 
that the duals of blocked edges make no cycles, i.e. form a forest (see figure 12 for an 
example) with no further restriction. The problem therefore amounts to enumerating 
forests on planar maps with tetravalent vertices, with a weight g per vertex and, say y 
per edge of tree in the forest. 

Imposing that both direction be blocked simultaneously is achieved in the mobile 
language by requiring that every (bivalent) black vertex on the mobile has either or 2 
incident marked edges. This constraint results in a slight modification of eq.(3J.), now 
replaced by: 

Qo{z) = -+[Q.{z)]^ + yQ.{z) 



Q.(z)=z + g [Qliz)] 
Q.{z)=gQl{z) 



R = l + g 



R 
J 



which amounts to a splitting of the operator Q,{z) of eq.( |3.1|) into Q»{z) and Q»{z) 

respectively corresponding to sub-mobiles with an un-marked and a marked root edge. 
The first three equations now reduce to 

Qoiz) = j + z + gyQoizf (3.9) 



from which we deduce the equivalent of eq.( p.7|) for the generating function R: 

n'.niin -f 1 ! 

n>l ^ ^ 

For instance, the first few terms in the expansion of R as power series of g read 

R = l+3g+6g^{3 + 5y) + 15g^{9 + 30y+28y'^) + 18g'^{63+315y 

+ 570|/2-f385|/^) + 126/(81 + 540y-t-14402/2 + 1855y^ + 100l/) (3.11) 
-f36£?^(2673 + 222752/ + 78300y2 + l469702/^ + 149884y4+68068y^)+C(/) , 

from which we read off the generating functions for forests on tetravalent maps (with 
distinguished face and unoccupied edge) with up to 6 vertices. In appendix A, we show 



20 



how to reinterpret eq. ( p.lO| ) in the language of even-valent maps without matter. Note 
that equation ( p.lO|) is equivalent to eq. (|3.7|) upon substituting g — > g{l + y)^. The 
origin of this extra factor is also explained combinatorially in appendix A. 

Let us now investigate the analytical behavior of as a function of g for a fixed 
finite positive value of y (note that each term of the series expansion of R in powers of 
g is a polynomial in y). Introducing the function 

n\n\{n + 1)! 

n>l ^ 

we can rewrite eq.( p.lO|) in a parametric form as 

u{y- F{u)) 

9 = 2 

J (3.13) 

R 



y - F{u) 



where we have introduced for convenience a parameter u = gyR. In this parametriza- 
tion, u must lie between and 1/27, which is the radius of convergence of F{u). How- 
ever, the first equation determines implicitly u in terms of g only in a range in which 
dg/du does not vanish, namely u G [0,w*) where is the smallest positive solution of 
dg/du{u*) = 0. This last condition may be rewritten as 

y = F{u*) + u*F'{u*) (3.14) 

which has a unique solution < u* < 1/27 for any finite positive value of y since 
the function F{u) + uF'{u) is strictly increasing and maps [0, 1/27) to [0, +oo). This 
implies that u{g) has a finite radius of convergence g* = u*{y — F{u*))/y'^. Note that 
d'^g / dv?{u*) < which implies that the singularity of u{g) has of the form {u* — u{g)) ~ 
{g* — g)^^^- Finally, we easily see that the function R{u) in ( |3.13| ) is regular in the range 
[0,w*] (as F{u) < F{u*) < y in this range) and dR/du{u*) > 0. Hence R{g) has the 
same radius of convergence g* and the same square root singularity. This corresponds 
to a critical exponent 7 = —1/2 (defined by R{g*) —R{g) ~ {g* — g)~^) characteristic of 
the universality class of so-called pure gravity, which is also found for quadrangulations 
without blocked edges. 

Finally, a different singular behavior is obtained by taking the y 00 limit in eq. 
( p.lOP , keeping a = gy fixed. In this limit, we readily see that R = 1 + {l/y)F{a) + 
0(l/?/2) where F{a) can now be interpreted as the generating function of maps with 
tetravalent vertices and equipped with a spanning tree, i.e. a single connected tree 
passing via all the vertices of the map, and with a weight a per vertex. The map has as 
before a distinguished face (origin face) and a distinguished unblocked edges separating 
two faces at distance m and m + 1 from the origin. We simply interpret the coefficient 
of a"" in F{a) as 

(3n)! (3n)! (2n + 2)! , ,1 , , 

^ X7 — ,,,, X (n + 2) X - (3.15) 



nlnl{n + iy. n!(2n + l)! {n + l)l{n + 2)l 
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Fig. 13: A rooted ternary tree (left). Upon matching leaves by pairs via 
a set of non-crossing arches, we obtain a tetravalent map endowed with a 
spanning tree, and with a distinguished oriented edge. 

expressing, as displayed in figure 13, that our tetravalent map with a spanning tree is 
made of a rooted ternary tree with n inner vertices and 2n + 2 leaves (first factor), 
completed with a system of n + 1 arches connecting these leaves (second factor) and 
forming the unvisited edges, with a choice of one among n + 2 faces (third factor). The 
last factor 1/2 is because the rooting of the ternary tree amounts to distinguish an 
oriented unvisited edge (that emerging from the root) while, in F, we distinguish only 
an un-oriented edge (here all edges are of type m m + 1, irrespectively of the choice 
of origin). For a ^ a* = 1/27, we now have the singular behavior F{a*) — F{a) ~ 
[a* —a)Log{a* — a), hence a critical exponent 7 = —1. The same exponent was obtained 
in Ref.|jl^ for a slightly different model. 

The above analysis applies straightforwardly to the case of section 3.1 by substitut- 
ing g g{l + y)'^. Among other possible direct applications of eq.( p.l4|) , one may also 
consider the case of arbitrary Eulerian maps with a maximal number of blocked edges. 
This case is discussed in detail in appendix B. 

4. Applications to maps with hard particles 

4.1. Particles with exclusion rules on maps 

Beside the direct applications already mentioned, the main interest of Eulerian 
maps with blocked edges is that they give access to the statistics of a number of quite 
involved "matter" models of maps equipped with interacting particles or spins. More 
precisely, we shall focus on the quite general case of Eulerian maps with particles sat- 
isfying a nearest neighbor p- exclusion rule. This rule is expressed as follows: each face 
of the map may be occupied by a number of identical particles with the restriction that 
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there be at most a total of p particles on any pair of adjacent faces WM. For p = I, the 



exclusion rule is nothing but the celebrated hard-particle condition that two particles 
cannot occupy identical nor adjacent sites. In general, in addition to the usual weights 
Qk and Qk per white and black /c-valent face, we attach the weight Zi to each face oc- 
cupied by i particles (with zq = 1). All these models are known to undergo various 
crystallization transitions to ordered phases in which the density of particles on black 
and on white faces are different (see Ref.[|15| for details). 




(a) (b) 

Fig. 14: A piece of map (a) with blocked edges and with some maximal 
domain (inside the dashed line) which is not attainable from the origin of 
the map, supposed to lie outside of the figure. All potentially incoming 
edges are blocked so that no path can enter the domain. In the dual 
representation (b), the boundary of the domain is a loop of length 2L of 
alternating black and white vertices, each carrying a number of particles, 
z = 0, . . . , 2L — 1 (and with n2L = ^o) and with every other edge blocked. 
We have n2i-i + n2i > p + 1 for i = 1 . . . , L, hence there exists some 
io G {0, . . . , L — 1} for which n2io + ?^2io+i > P + 1- The corresponding 
edge may be freely blocked or not, thus contributing a weight 1 + y — 0. 



The p-exclusion rule may be implemented as follows: we first start by relaxing it 
and demand only that there be at most p particles on each face. We then introduce on 
top of any such particle configuration a set of blocked edges on the map, with a weight 
y per blocked edge, and with the two following restrictions: 

(i) as before, we demand that blocked edges satisfy the global connectivity constraint; 

(ii) we furthermore restrict the configurations of blocked edges by allowing only for the 
blocking of those edges which violate the p-exclusion rule for the particle configu- 
ration at hand. 

As we already know, restriction (i) is crucial to have a mobile description of the config- 
urations at hand. Restriction (ii) is a new constraint but it will be completely transpar- 
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ent in the mobile formalism. Finally, the p-exclusion rule is restored by simply taking 
y = —1 in the various generating functions. 

Indeed, note first that whenever the configuration of particles satisfies the p- 
exclusion rule, we cannot add any blocked edge and it receives the correct weight (ir- 
respectively of the value of y). If the configuration of particles violates the p-exclusion 
rule, let us show that the contributions of all blocked edge configurations allowed on the 
map add up to 0. Clearly, if we relax the global connectivity constraint (i), any violated 
edge can be blocked or not, thus receiving a weight 1 + y = so that the contribution 
of all blocked edge configurations vanishes trivially. To prove that the contribution of 
those configurations that satisfy (i) also vanishes, we therefore need only prove the van- 
ishing of the complementary contribution of those blocked configurations that do not 
satisfy (i). To this end, assume the existence of a maximal connected domain which is 
not attainable from the origin of the map and consider the boundary of the dual of this 
domain (see figure 14- (a)). It consists of a loop of alternating black and white vertices 
connected via a succession of edges dual to alternatively in- and outgoing edges of the 
original Eulerian map (see figure 14- (b)). Clearly, all in-going edges must be blocked to 
prevent entering the domain. We simply have to show that at least one of the outgoing 
edges violates the p-exclusion rule so that we may freely decide to block it or not, thus 
creating a weight 1 + y — causing the whole contribution to vanish. As all in-going 
edges are blocked, they violate the exclusion rule, thus the total number of particles 
around the loop is larger or equal to {p + 1)L if the loop has length 2L. Assuming 
by contradiction that none of the out-going edges violates the exclusion, then the total 
number of particles around the loop would not exceed pL. This forms a contradiction. 

For illustration, we have displayed in figure 15 the allowed blocked edge configura- 
tions satisfying (i) and (ii) above for p = 1 and for a particular configuration of particles. 
Their total contribution vanishes when y = —1 as the particle configuration at hand 
violates the 1-particle exclusion rule. 

We conclude that the nearest-neighbor interaction of the particle model with p- 
exclusion rule can be completely undone at the expense of introducing blocked edges 
on the map. This will allow for a straightforward derivation of recursive equations 
determining the generating function of the model. 

4-2. Generating functions 

Let us consider the generating function G of Eulerian maps with particles subject 
to the p-exclusion rule, with a distinguished origin vertex (pointed maps) and, for con- 
venience, with also a distinguished edge. Thanks to the bijection of section 2, G may 
be expressed in terms of various generating functions for rooted sub-mobiles. More pre- 
cisely, the mobiles that we have to consider now are unrestricted mobiles supplemented 
by particles living on their black and white unlabeled vertices, and with at most p par- 
ticles per vertex. Given this particle decoration, the only new constraint inherited form 
the condition (ii) of previous section is that a fiagged edge may be marked only if the 
total number of particles on its two endpoints exceeds p strictly. By a direct general- 
ization of our notations of section 2.4, we may now introduce the following generating 
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Fig. 15: A sample configuration of pointed Eulerian map with 4 particles 
(top left) violating the p-exclusion rule, here with p = 1. We have listed 
all possible blocked edge configurations satisfying the global connectivity 
constraint (i) and the condition (ii) that we may block only edges that 
violate the 1-exclusion rule. Weighting each blocked edge hj y — —1, the 
contributions of the 12 configurations add up to 0, as it should. 

functions: 

- R denotes as in eq. ( |2.9| ) the generating function for sub- mobiles rooted at a labeled 
corner, and with arbitrary (but fixed) root label; 

- B^^^ (resp. B^^^) denotes the generating function for sub-mobiles rooted at some un- 
marked (resp. marked) fiagged edge incident to a black vertex carrying i particles 
and with arbitrary (but fixed) labels on the left and on the right of the edge differing 
by i (with £ > in the un- marked case); 

- finally, W^j (resp W^^) denotes to the generating function for sub- mobiles rooted 
at some un-marked (resp. marked) flagged edge incident to a white vertex carrying 
i particles and with arbitrary (but fixed) labels on the right and on the left of the 
edge differing by i (with again £ > in the un-marked case). 

By use of a straightforward generalization of eqs.( p.ld|) and (|2.11|) with y = —1, the 
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generating function G is now given by: 

i=o j=o e>o i=o j=p+i-i e 

Collecting the generating functions of the different sub-mobiles that can be attached to 
a white (resp. black) vertex carrying i particles, we are led to introduce p + 1 Laurent 

series Q^o\z) (resp. Qi^\z)), i = 0, . . . ,p, defined as: 

R ' " 



j=0 £>0 j=p+l-i I 



(4.2) 



j=0 £<0 j=p+l-i I 



All generating functions are now fixed by the recursive relations: 

i? = Y^, L = 5^^.X^(g«(z))'"\. 

k i=0 

Bf=z, $^^fc(Qi^^(^))'"\^, ^>0 

k 

Bf = -z. $^^.(Qi^^(^))'"\. (4.3) 

k 



Wf^ = -z. J29k {Q'i\z) 



k-i 



In particular, note that ^ = -Sf ^ and W^J:} = -W^^] for £ > 0, so that eq. ( P1 ) 
simplifies into 

« = -R+EEE<»'5-E E E^f'w'ii' (4.4) 

i=0 j=0 ^>0 i=0 j=p+l-i e<o 

while we may write eqs.( [4.2| ) as 

p—i 



^"w = 7+EE<'~'+ E E^f-' 

j=0 £>0 j=p+l-i£<o 



j=0 £<0 i=p+l-i ^>0 



26 



(4.5) 



i+j < 2p 

iO |) 



i+j ^ p 

io-^ — mj 



i+j ^2p 




p+1^ i+j 4 



2p 




^^^^^^^^ • 




'J 



no restriction on labels 



m<n 



(a) 



(b) 



Fig. 16: Simplification of tfie crossing rule wfien y = —1. Here i and j 
denote tlie numbers of particles occupying the white and black vertices. 
The original crossing rule (a) of section 2.2 displayed in figures 4 and 5 can 
be equivalently replaced by a new crossing rule (b) in which the allowed 
variation of labels m, n is strictly correlated to the satisfaction (m > n) or 
violation (m < n) of the p-exclusion rule. 



In particular, we need only consider generating functions for i < (and VF^ for 
i > 0). This is a direct manifestation of an obvious simplification of the crossing rule 
for labels on mobiles when y = — 1, as displayed in figure 16. Any fiagged edge whose 
endpoints have a total number of particles larger or equal to p + 1 and with increasing 
label clockwise around the white endpoint can be marked or not on the mobile, resulting 
in a weight 1 + y = 0. The mobiles that survive in practice are therefore mobiles whose 
fiagged edges are of the following two types only (see figure 16-(b)): 

- unmarked fiagged edges whose endpoints have a total number of particles smaller 
or equal to p and with increasing label clockwise around the white endpoint; 

- marked flagged edges whose endpoints have a total number of particles larger or 
equal to p+1 and with strictly decreasing label clockwise around the white endpoint. 
These edges furthermore receive a weight —1. 

As in section 2.4, generalizing ( |2.14| ), we may summarize the equations above in a closed 



27 



system: 



p—i 



k j=0 
p—i 



k j=0 

p 



k j=i 



k-l 



k-l 



k-l 



V 



~l~ k j=p+l — : 

P 

k j=p+l— i 



R 



Q^^\z] 



k-l 



k-l 



<0 



>0 



(4.6) 

where the notations [■]<o and [•]>o stand respectively for the extraction of the negative 
and positive power parts in the Laurent series of z. 



4.3. Detailed study of Eulerian triangulations with hard particles 

As an example, let us study in detail the case of Eulerian triangulations with hard- 
particles, i.e particles subject to the 1-exclusion rule. This case corresponds to having 
9k = 9k = 9Sk,3 (as there are as many black and white triangles in an Eulerian triangu- 
lations, there is no need to introduce different weight factors for these faces). In a dual 
formulation, this model is equivalent to that of so-called bicubic maps with hard parti- 
cles studied in Ref. ||15|| by random matrix techniques and in Ref. ||14|| by combinatorial 
techniques based on blossom trees. We refer to these references for a description of the 
physical properties of the model. By a simple look on eqs. ( |4.5D and ( ^.tj| ), we imme- 
diately deduce that the only non vanishing generating functions in this case are 



B^P for I = 
Vrls^ W^f ^ for 



2, 5, 8 and sj^^ for 



—2, —5, —8 and wi^^ for £ = 1,4, so that we may write eq. 



-1, —4 together with their "white" counterparts 

) as: 



Q^^\z) = ^ + Bi'^z' + B^'^z' + Bi'^z' + Bi'^^z' 

D r(1) r(1) 

Z^ Z'^ z^ z^ 

Qi^^(.) = . + ^^ + W^W. + T^iV. 



Equations (^4.3|) now translate into: 



= 9 = gR- 

Bi'^ = ^.1 ((l + Wi'^y + 2W^_'M'^^ W^:^ =9z^^{R + 5^)' + 23^'^^) 

= -2gz,W^_'^ (l + Wi'^) Wi'^ = -2gz,B'i^ {r + B^-1) 

B^y, = -gz,{w^^^y wi'^ = -gz,{Bi^^y 

(4.8) 

as well as expressions for 55"*^^ B^\ ^^-5* and W^g that we shall not use. This allows 
to express these generating functions in terms of R only, namely: 

Bi'^=g 



(1) ^ gzijl - 2g^z,R^ - Sg^zjR^ - Sg^zfR'') 
2 {l + 2g^ziR)^ 



(4.9) 



l + 2g^ziR 



B^ll = -g^z.R' , 

while W^^l = R^Bf\ W^_}1 = R^B^^^\ W^^^ = B^}l/R and w[^^ = B^_}l/R\ Finally, 
the generating function R itself is given by the recursion relation: 



R = l+gR{2R (^Bi'^ + S^^^) + 2z, (i? + B^'l) B^^^ 



(4.10) 



(1 + 2£?2;2li?)2 

In particular, expanding R in powers of g leads to: 

R = l + 2g^{l + 2zi)+4:g^{2 + 8zi+5zf)+4:g^{10+60zi+89zl + 28zf) + 32g^{7+56zi 
+ 135zj + 107zf+21zf) + 16g^^{84: + 840zi + 2828zl+3808zf + 1911zf + 264:zl) 
+ 64t/^2 (^1^2 + 1584;2i +6870^? + 13320^^ + 1 1629^^ +4088z^ +429z^ ) + 0{g^^) . 

(4.11) 

The desired generating function G for maps can be expressed in terms of R as: 
G = R + Bi'^W^'^ + + Bi'W^^ - B^_'lwi'^ - B^^M^ 

{l + 4g^ziR^ +g^{R + 6ziR)-g'^zfR^ -20g^zfR'' -20g^^zfR^) (4-12) 
~^ {l + 2g^ziR)^ ■ 

In particular, expanding G in powers of g leads to: 

G = l + 3g^{l+2zi) + 12g^{l+4:Zi + 2zl) + 15g^{4: + 24:Zi + 33zf + 8zf)+48g^{7+56zi 
+ 1302^ + 922^ + 144) + 168£^^°(12 + 1202i+3952f + 500;zf + 2204+242f) 
+ 144^^2 (88 + 105621+45122^+84162^ + 68012^+2080^5^ + 1762^) +C(£?^^) , 

(4.13) 
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which displays the generating functions of Eulerian triangulations with hard particles 
and with 2, 4, . . . , 12 triangles. Note that here, the triangulations have a distinguished 
origin vertex and a distinguished edge. It is more usual to consider so-called rooted maps 
with a distinguished edge only. As there are n + 2 vertices in an Eulerian triangulation 
with 2n faces, this is done here by dividing out the coefficient of g^^ by an overall factor 



Fig. 17: An example (a) of quadrangulation with Ising up/down spins 
on its faces and the equivalent Eulerian map with bi- and tetra-valent 
black or white faces (b) . The original spin up (resp. down) faces translate 
into empty black (resp. white) tetravalent faces, while a bivalent face is 
introduced between any pair of adjacent faces with equal spins to restore 
the bicolorability. The same configurations are alternatively selected from 
all possible Eulerian maps with bi- and tetra-valent faces by adding a 
particle on each bivalent face and imposing the 1-particle exclusion. This 
indeed forces the bivalent faces to be isolated and to separate tetravalent 
faces of identical color. 

As a final application, let us consider the case of quadr angulations with Ising spins, 
i.e maps with tetravalent faces only and with spins up or down living on these faces. The 
Ising coupling is introduced by a weight, say zi per edge separating faces with the same 
value of the spin, while edges separating faces with opposite spins receive the weight 1 
instead. The model may be realized in terms of Eulerian maps with hard-particles as 
follows (see figure 17): consider Eulerian maps made of: 
- black or white tetravalent faces that carry no particle; 



(n + 2). 



4.4. Ising model revisited 




(a) 
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- black or white bivalent faces that, on the contrary, are forced to carry a particle, 
weighted by zi. 

Assume now that the particles are subject to the 1-particle exclusion rule. This imposes 
that two bivalent faces cannot be adjacent to one another so that any bivalent face on 
the map necessarily lies between two tetravalent faces. If the bivalent face is black, 
the two neighboring tetravalent faces are necessarily white (and conversely) so that, 
by squeezing all bivalent faces into single edges, we obtain a quadrangulation with 
black and white faces and with an effective weight factor zi per edge separating two 
faces of the same color. This is nothing but the desired Ising model by interpreting 
the black (resp. white) faces as carrying spins up (resp. down). The Ising model 
thus corresponds to considering Eulerian maps with particles subject to the 1-particle 
exclusion rule, with Qk = dk = gSk,^ + zidk,2 and with the additional restriction that 
bivalent faces must carry particles while tetravalent faces must be empty. Again this 
last restriction is transparent in the mobile formalism as it simply means that the black 
and white vertices of the mobile must carry a particle or not according to their valence 
on the mobile. By a simple look on eqs. (|4.5|) and ([4.3| ), we immediately deduce that 



now, the only non vanishing generating functions are ^ and ^ for £ = 1,3 and ^ 
for i = —1, —3, together with their white counterparts W^^^ and W^^'^ for i = —1, —3, 



and Wf^'^ for £ = 1,3. We may thus write eq. (|4.5|) as: 



p r(i) r(i) 

Q^}\z) = I + Bfh + Bfh' + ^ + ^ 

(0). . W^'l W^'l W^'l W^'l 

Z Z'^ z z-^ 

z z-^ 



(4.14) 



Equations (|4.3| ) now translate into: 





= 9 




Bf 


= 3g {W^^l + W['^) 




^(1) 


= -z^W^'i 




B^\ 


= -Z^W'^'I 






= z^W^'^ 






= .1 (i + wi'') 





39 R^ (Sf^+Sf) 



-z^B['\ 



(4.15) 

zi (r + b^'I) 
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from which we deduce ah generating functions in terms of R only: 

(0) _ 3gziR 

^1 ~ 1 Q„^1 



l-3g{l-zl)R 
^(1) zi-3gziR (4.16) 



1 -l-3g{l-zf)R 
^ll=-9z,R' 

(1) 'igzlR'' 



TD\ 

1 - 3^7(1 -^?)^ 

and = R^Bg^^ for aU the values of £ and i mentioned above. As for R itself, it now 
satisfies the recursion relation: 



R=l + R (^3gR' (Sf + B^'^) + S^/i? + S^^)' + z^B['^ 

^^2zj - 6g^R^{l - z\f + %'^R^{\ - zlf + g R{1 - + 3zt] 



(4.17) 



Note that this equation may be rewritten as 

— % = 1 + 3g'{l - zl)R' + , , (4.18) 



(l-zi) ^ ^ \-zi{\-3g{\-z{)R)'' 



which is precisely the form found in Refs. |T^], [0 and |TJ]. The first terms of the 
expansion of R in powers of g read: 

R = l+6£?z?+3/(l+8z? + 15zf) + 18£?^2^(ll + 28z? + 21zf)+27/(l+3l2? 

+ 229;2t + 285z^ + 126z?) + 54/2f(97+907zf + 252l2t + 19292^+594;2^) (4.19) 
+8l£^^(4 + 2792?+4833z^ + 199582^+30678z? + 164192}° + 386l2}2)+O(/) . 

The generating function G for maps is now expressed in terms of R as: 

G = R + Bfw'-^l + Bfw^^l + SJ^VI^) + S^^VS,) + Bfw^ll + 
~ B^^XW"^^ - b'^Hw^^^ 

_ \^\^g'^R^{\-3zl)-^gR{\~2zl)^%''R''{\-zlf{\-3zl)-^g''R^{\-^zl^3z\) 

{\-3g{\- z^YRY 

(4.20) 
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and expanding G = G{g, zi) in powers of g leads to: 

G{g,zi) = l + 12gzf+4:g^(l + 12zl + 18zf) + 180g^zl{2 + 5zl+3zf) + 18g'^{2+85zl 

+ 624:zj + 693zf+252zf) + 756g^zf{12 + 119zl + 312zj + 207zf + 54:zf) 

+432/ (1 + 91^^ + 1642zj + 6681zf+9450zf+A356zl^+891zl^)+O{g^) . 

(4.21) 

The coefficient of g'^ in G is the generating function of maps with a distinguished 
origin vertex (among n + 2) and with a distinguished edge in the original Eulerian map 
formulation with bivalent faces. This means that on the associated quadrangulations, 
edges separating two faces of the same color can be marked twice. The number of such 
edges is given by the power of ^i, so that, in the end, coefficient of g^z\^ in ( ^4.21| ) receives 
a weight (n + 2)(2n + 2p) for all possible markings rather than the usual factor 4n for 
rooted quadrangulations with a distinguished oriented edge only. We can easily restore 
this factor and obtain the generating function H = H{g, zi) for rooted quadrangulations 
with Ising spins, with expansion: 

H{g, zi) = 4g^ C da C d(3^ {Gi(3a^g, az,) - l) 
dgJo Jo 

^4:gzl+2g^{l + 8zf+9zt) + 108g^zf{l + 2zl+zt) + 12g^{l+34:zf 

+ 208zi + 198zf+63zf) + 216g^zl{10 + 85zl + 195zi + 115zf+27zf) 

+ 54:g^{2 + 156zf + 2463zj+8908zf + 11340zf+A752zl^ + 891zl^)+O{g^) . 

(4.22) 



5. Conclusion 

In this paper, we have considered the combinatorics of the general class of Eule- 
rian maps with blocked edges. This includes a number of matter statistical models on 
(Eulerian or unrestricted) maps whose local interactions are mediated by the choice of 
the blocked edge configurations. All these models are amenable to a coding by mobiles, 
and we obtained the general recursive system of equations that determines the associ- 
ated generating functions. We gave a number of explicit direct or indirect applications. 
In particular, this provides a new combinatorial approach to the problem of mutually 
excluding particles, a set of models that gives access to most relevant RCFT critical 
points. 

The general structure of the recursive equations above involve Laurent series Q»{z) 
and Qo{z) or, more generally, operators Q, and Qo which appear as transfer opera- 
tors describing the evolution of labels around vertices of the mobile. Remarkably, these 
operator bare a striking similarity with the "position" operators standardly introduced 
in the chain-interacting multimatrix models solved by bi-orthogonal polynomial tech- 
niques. In this framework these linear operators correspond to the multiplication by an 
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eigenvalue of some random matrix, usually expressed in the basis of orthogonal poly- 
nomials. This similarity was already observed in the simpler mobile or blossom tree 
descriptions of maps without matter and seems to indicate a deep connection between 
the existence of bijections with trees and the solvability of the corresponding matrix 
models via orthogonal polynomials. This is yet to be understood. 

At a more practical level, the recursive equations that we face involve in general 
truncations of the Laurent series Q,{z) and Qo{z) to both their positive and negative 
parts. In the applications that we discussed, we have explored only particular cases 
in which either we may get rid of the truncations (section 3) or only a finite number 
of coefficients in the Laurent series survive (section 4). In the first case, this leads to 
simple and easily solved algebraic equations for Q,{z) and Qo{z). In the second case, 
we end up with a closed algebraic system for the finitely many surviving generating 
functions. Note that we have not explored the most general case where infinitely many 
coefficients would survive, as we do not really know how to solve the system of equations 
in this case. Finding the corresponding solution could give access to more general (non 
necessarily algebraic) critical behaviors. Such behaviors have been observed for other 
types of statistical models such as loop models [0, Potts models [|T^ or the 6- vertex 
model p^. The latter correspond to matrix models that were solved by loop equation 



or saddle point techniques, as no orthogonal polynomial techniques were available. We 
may wonder whether our general class of Eulerian maps with blocked edges could deal 
with some of these problems. 

To conclude, the general framework of our approach is a modification of the natural 
notion of geodesic distance on the maps, correlated with the presence of matter. This 
is a quite natural principle of general relativity, here transposed at a discrete level. In 
the present construction, the distance was modified by introducing blocked edges but 
we may imagine other mechanisms. Note finally that our formalism allows in principle 
to keep track of this modified distance (by keeping track of the labels in generating 
functions), thus giving access to refined intrinsic geometrical properties of maps with 
matter. 
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Appendix A. Combinatorial interpretation of eqs.( |3.7D and (|3.10|) : forests 
on tetravalent maps as decorated even-valent maps 



Equations (|3.71 ) and (|3.10|) were obtained by use of our mobile formalism but can be 



given a simple interpretation. These equations are strongly reminiscent of that obeyed 
by the generating function (still denoted by R) of planar maps with a distinguished face 
and a distinguished edge and with vertices of arbitrary even valence 2k, counted with 
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weight V2k (see 0): 

i?=i+x:^2.(\^)i?^ (A.l) 

Eq.( p.lO| ) is recovered by setting 

k-i,,k-2 (3/c-3)! 

for all /c > 2 and f2 = 0. This relation may be understood as follows: any connected 
component of the forest with, say k — 1 tetravalent inner vertices (hence k — 2 inner 
edges and 2k leaves) can be contracted into a single 2/c-valent vertex by squeezing all its 
inner edges. Conversely, any 2/c-valent vertex can be expanded into a ternary connected 
tree with k — 2 edges in exactly {3k — 3)!/((/c — l)!(2/c — 1)!) ways. The enumeration of 
forests on maps with tetravalent vertices is therefore equivalent to that of maps without 
forests but with vertices of arbitrary even valences, with the weight factor V2k above 
simply accounting for the appropriate degeneracy factor in the squeezing process. 

Concerning eq.( |3.7| ), it differs from eq.( |3.1Cl| ) only by an extra factor (1 + y)^"^ which 
can be simply understood from a 1 to 2^" correspondence, for pointed quadrangulations 
with 2?! edges, between allowed configurations of blocked edges where the blockings are 
necessarily in both directions and allowed configurations of blocked edges where block- 
ings in a single direction are also allowed. The correspondence is as follows: starting 
from a configuration with edges that are either not blocked or blocked in both direc- 
tions, we consider for each edge the orientation that leads from its endpoint with a larger 
distance to the origin to that with a smaller distance (on a quadrangulation, distances 
between neighbors have different parities so one must be strictly larger than the other) . 
For this particular orientation, we keep the blocked or unblocked nature of the original 
edge but for the reverse orientation, we decide arbitrarily to block it or not. This leads 
to 2^" configurations with the same distance on the map as the presence or absence of 
these last blockings clearly does not affect this distance. Collecting the y factors for 
these 2^"^ configurations clearly produces a term (1 -(- y)^". The reverse construction is 
as follows: starting from a configuration with possible blockings in a single direction, we 
first suppress all blockings oriented from a vertex at smaller distance from the origin to 
a vertex at larger distance, which leaves us with edges which are either not blocked or 
blocked in one direction only (one-way edges). We then transform these one-way edges 
into edges blocked in both directions. 



Appendix B. Maximally blocked Eulerian maps 

The contribution of Eulerian maps with a maximal number of blocked edges is 



obtained by taking in eqs.( 2.14 ) the limit y - 
finite for all k. In this case, we may write R 



00 keeping ak = Qky and ak = Qky 

1 + {l/y)Z + 0{l/y^) where Z is the 
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Fig. 18: A bicolored tree with as many blaclc leaves as wiiite leaves (left). 
There is a unique matching between black and white leaves that can be re- 
alized with non-crossing arches and is such that every face can be reached 
from the external face upon respecting the canonical orientation, i.e. cross- 
ing arches with their black endpoint on the right. 



generating function of maximally blocked Eulerian maps. From eq.( |2.14| ), we have the 
relations 



Z = y^ak 

k 

Qo{z) = - + Y,akQ'.-\ 

^ k 

Q.iz) = z + J2^kQ'.~'i< 



(B.l) 



These equations are recognized as the recursion relations for the generating functions 
of rooted planar trees whose nodes (inner vertices) have alternating black or white 
color, counted with weights per /c-valent white node and ctfc per /c-valent black node. 
More precisely, interpreting the parameter 2 as a weight per leaf adjacent to a black 
node (hereafter referred to as black leaf), while a weight 1/z is attached to each leaf 
adjacent to a white node (white leaf), we see that {Q,{z)/z) — 1 (resp. zQo{z) — 1) 
is the generating function for such trees with a distinguished white (resp. black) leaf. 
Rewriting the first line of eq.( |B.l| ) as 



'Q.{z) 



1 



(B.2) 
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we see that Z can be interpreted as the generating function for bicolored trees with a 
distinguished white leaf and constrained to have as many white leaves as black leaves. 
This should not be a surprise. Indeed, in a maximally blocked Eulerian map, the set of 
edges dual to blocked edges forms a spanning tree of the dual map, which is naturally 
bicolored. Each non- blocked edge on the dual map connects nodes of opposite colors: 
cutting it into two half-edges terminating at leaves results into a bicolored tree with as 
many white as black leaves. The original map has a distinguished non-blocked edge, 
which turns into a distinguish white leaf. Conversely, given a bicolored tree with as 
many white leaves as black leaves, there is a unique way to connect each black leaf 
to a white one via non-intersecting arches so that each face of the resulting map can 
be reached from, say the external face upon respecting the canonical orientation, i.e. 
crossing edges only with the white endpoint on the left and the black endpoint on the 
right. This can be done by an iterative matching algorithm: we follow the contour of 
the tree clockwise and match each black leaf immediately followed by a white leaf to 
this latter. We then iterate the process by ignoring the already matched leaves until all 
leaves are matched (see figure 18). The resulting map is clearly dual to a maximally 
blocked Eulerian map, with a distinguished origin vertex (dual to the external face) and 
a distinguished non-blocked edge of type m m + 1 (dual to the edge connecting the 
distinguished white leaf), and whose blockings satisfy the global connectivity constraint. 
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